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Abstract 

In  this  paper  we  present  a new  formulation  of  a large  class  of  phase-held  models,  which  describe 
solidification  of  a pure  material  and  allow  for  both  surface  energy  and  interface  kinetic  anisotropy,  in 
terms  of  the  Hoffman-Cahn  ^-vector.  The  ^-vector  has  previously  been  used  in  the  context  of  sharp 
interface  models,  where  it  provides  an  elegant  tool  for  the  representation  and  analysis  of  interfaces  with 
anisotropic  surface  energy.  We  show  that  the  usual  gradient- energy  formulations  of  anisotropic  phase- 
field  models  are  expressed  in  a natural  way  in  terms  of  the  ^-vector  when  appropriately  interpreted. 

We  use  this  new  formulation  of  the  phase-field  equations  to  provide  a concise  derivation  of  the  Gibbs- 
Thomson- Herring  equation  in  the  sharp-interface  limit  in  three  dimensions. 

"Address  from  Jan.  1st  ’94:  Faculty  of  Mathematical  Studies,  University  of  Southampton,  Highfield,  Southampton,  S09 
5NH,  UK. 
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1.  Introduction 

The  presence  of  surface  energy  anisotropy  in  the  solidification  of  a pure  material  is  a phenomenon  of 
both  practical  and  theoretical  importance.  In  practice  it  leads  to  non-spherical  equilibrium  shapes, 
which,  if  the  anisotropy  is  sufficiently  pronounced,  may  exhibit  missing  orientations  (e.g.,  edges  or 
corners)  or  facets  (flat  surfaces).  Its  importance  has  been  recently  highlighted  by  the  suggestion  that 
it  may  be  responsible  for  the  selection  of  the  dendrite  tip  operating  state  in  the  growth  of  a dendrite, 
see  [1]  for  a review.  Computations  of  dendritic  growth  using  the  phase-field  method  in  two  dimensions 
also  suggest  that  surface  energy  anisotropy  is  an  important  effect  in  dendritic  growth  [2,  3,  4].  Various 
ad  hoc  modifications  to  phase-field  models  have  been  proposed  to  represent  surface  energy  anisotropy. 
However,  only  for  the  model  proposed  by  Kobayashi  [3]  is  the  corresponding  boundary  condition  for  the 
interfacial  temperature  in  the  sharp  interface  limit  known  and  then  only  for  an  interface  represented  as 
a one- dimensional  curve,  see  [5].  In  reality,  dendrites  are  three-dimensional  and  a sharp  interface  model 
requires  the  representation  of  their  interfaces  as  two-dimensional  surfaces.  A phase-field  model  can 
also  provide  a setting  for  computation  of  three-dimensional  dendrites,  see  [6].  In  order  to  assess  such 
computations  it  is  important  to  know  the  free-boundary  problem  approached  in  the  sharp  interface 
limit.  It  is  the  aim  of  this  paper  to  address  this  issue. 

In  this  paper  we  provide  a new  formulation  of  a class  of  phase-field  models  that  allow  anisotropic 
surface  energy  and  interface  kinetics  in  terms  of  the  Hoffman-Cahn  ^-vector  [7,  8].  We  exploit  this 
representation  to  show  that  in  the  sharp-interface  limit  these  phase-field  models  yield,  at  leading 
order,  a free  boundary  problem  in  three  dimensions  which  correctly  accounts  for  the  surface  energy 
and  kinetic  anisotropy.  In  a previous  paper  [5]  we  derived  the  corresponding  result  in  two  dimensions 
for  a specific  phase-field  model  using  a rather  lengthy  asymptotic  analysis.  In  contrast,  the  vector 
formulation  we  develop  here  allows  for  a concise  derivation  in  three  dimensions  of  a whole  class  of 
anisotropic  phase-field  models. 

Hoffman  and  Cahn  [7]  introduced  the  notion  of  a ^-vector  to  represent  the  anisotropic  surface 
energy  of  a sharp  interface  separating  two  phases.  It  provides  a particularly  elegant  and  useful  device 
to  describe  an  interface  with  an  anisotropic  surface  energy,  7(n)j  where  n is  the  unit  normal  vector  to 
the  interface.  It  yields  a succinct  form  of  the  equation  for  the  interfacial  temperature,  T/, 

Ti  = Tm-^Vs-(,  (1) 

which  can  be  shown  to  be  equivalent  to  the  classical  Gibbs-Thomson- Herring  equation;  here  we  have 
neglected  the  effects  of  interface  kinetics,  which  we  consider  in  more  detail  below.  Tm  is  the  melting 
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temperature,  L is  the  latent  heat  of  fusion  per  unit  volume  and  V5  is  the  surface  divergence  operator 
in  the  interfacial  surface  S applied  to  the  ^-vector.  The  surface  divergence  of  the  ^-vector  is  equivalent 
to  the  notion  of  generalised  mean  curvature,  as  discussed  by  Taylor  et  al.  [9,  10].  For  the  case  of  a 
two-dimensional  geometry  this  reduces  to  the  well-known  form  [11] 

T/  = 2M-^(7  + 7^^)iC,  (2) 

where  (p  is  the  angle  of  the  normal  n of  the  interface  to  a fixed  direction,  and  K.  is  the  local  curvature 
of  the  interface. 

Various  modifications  to  phase-field  models  to  allow  for  surface  tension  anisotropy  have  been 
suggested.  Danger  [12]  proposed  including  terms  of  squares  of  higher  derivatives  of  the  phase  field 
and  gave  an  example  leading  to  cubic  anisotropy.  Caginalp  and  Fife  ([13],  see  also  [14])  suggested 
that  the  square  gradient  term  that  appears  in  the  free-energy  functional  could  be  replaced  by  a more 
general  quadratic  form  with  different  coefficients  in  each  coordinate  direction.  A different  approach  is 
to  include  nearest  neighbour  interactions  in  a discrete  formulation  of  a diffuse  interface,  see  Cahn  and 
Kikuchi  [15].  More  recently  Kobayashi  [3]  and  Wheeler  et  al.  [2]  have  allowed  the  coefficient  of  the 
gradient  energy  to  depend  on  the  local  orientation  of  the  gradient  of  the  phase  field.  An  asymptotic 
analysis  corresponding  to  the  sharp  interface  limit  in  two  dimensions  was  conducted  by  McFadden  et 
al.  [5]  to  show  that  the  appropriate  form  of  the  boundary  condition  for  the  interfacial  temperature 
was  obtained. 

In  section  2 we  briefly  review  the  Cahn-Hoffman  theory  of  the  ^-vector  for  a sharp  interface.  In 
section  3 we  show  how  the  f- vector  arises  in  a natural  way  for  a class  of  anisotropic  phase-field  models. 
In  section  4 we  examine  a sharp  interface  limit  of  the  anisotropic  phase-field  equations  and  recover 
the  Gibbs-Thomson- Herring  equation.  Conclusions  are  given  in  section  5.  Some  technical  points  are 
relegated  to  appendices,  which  contain  summaries  of  the  relevant  diflferential  geometry  and  asymptotic 
expansions  that  are  employed  in  the  body  of  the  paper. 


2.  The  Hoffman- Cahn  vector 

Here  we  briefly  review  the  elements  of  the  theory  of  the  ^-vector  given  by  Hoffman  and  Cahn  [7]. 

For  an  isotropic  surface  energy  7,  the  equilibrium  shape  of  a solid  particle  surrounded  by  its  melt 
is  described  by  a form  of  the  Laplace- Young  equation 

ifC  = AF,  (3) 
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where  K is  the  local  mean  curvature  of  the  solid/liquid  interface  and  the  constant  AF  is  the  difference 
in  the  bulk  free  energies  between  the  liquid  and  solid  phases.  The  Laplace- Young  equation  can  be 
derived  by  minimizing  the  surface  energy  subject  to  a volume  constraint.  The  equation  may  also  be 
written  in  terms  of  the  surface  divergence, 

V5  • (7^)  = (4) 


where  n is  the  local  normal  to  the  interface.  Properties  of  the  surface  differential  operator  V 5 are 
reviewed  in  Appendix  A;  the  property  relevant  here  is  that  V5  • n = /C. 

In  the  anisotropic  case,  the  surface  tension  of  the  solid/liquid  interface  depends  on  the  local 
orientation  of  the  interface,  which  can  be  expressed  as  a functional  dependence  of  the  the  surface 
tension  on  the  local  unit  normal  to  the  interface,  viz.,  7 = 7(n).  Cahn  and  Hoffman  show  that  from 
the  function  7(71)  a vector  ^ can  be  defined  such  that  the  anisotropic  version  of  Eq.  (4)  which  governs 
equilibrium  shapes  has  the  form 

Vs-i  = AF-,  (5) 

the  ^-vector  reduces  to  the  vector  jn  if  7 is  constant.  Cahn  and  Hoffman  give  a formal  definition  of 
the  ^-vector  in  the  form 

( = V{n{e,(f));  (6) 


here  we  have  expressed  the  normal  vector  in  terms  of  the  spherical  coordinates  6 and  cp,  and  have 
written  = 'y{n{6,ip)).  The  gradient  appearing  in  Eq.  (6)  is  explicitly  given  by  the  expression 


z*  . ^ d'j  1 dj  _ 

t — j L0  

^ ^ de  sinddcp^' 


(7) 


where  r,  6,  and  (p  are  unit  vectors  in  the  coordinate  directions.  It  should  be  emphasised  that  the  unit 
vectors  f,  6,  and  (p  belong  to  a coordinate  system  associated  with  a specific  point  on  the  interface,  and 
are  used  to  describe  vectors  based  at  that  point.  In  particular,  in  this  coordinate  system  the  normal 
to  the  interface  is  given  by  n(0,  ip)  = </)),  so  that  the  normal  component  of  the  ^-vector  is  given  by 

7,  as  in  the  isotropic  case.  In  the  anisotropic  case,  ^ also  has  components  that  lie  in  the  plane  tangent 
to  the  surface,  which  is  spanned  by  the  unit  vectors  6 and  (p. 

In  the  above  definition  the  vector  is  parametrised  by  the  spherical  coordinates  of  the  normal 
vector  to  the  interface,  with  ^ If  the  equilibrium  shape  is  smooth  and  convex,  then  the 

angles  6 and  ip  can  themselves  be  used  in  a parametric  description  of  the  interface  itself;  that  is, 
the  interface  can  be  expressed  in  the  form  x = xijd^ip)^  where  6 and  ip  range  over  the  unit  sphere. 
Although  this  appears  to  be  a highly  implicit  description,  since  the  interface  shape  is  parametrised  in 
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terms  of  its  normal,  it  leads  to  an  explicit  expression  for  the  interface  shape.  This  follows  from  the 
identity  (see  Appendix  A) 

V5-x  = 2,  (8) 


stating  that  the  surface  divergence  of  the  position  vector  of  the  interface  is  two.  Comparing  Eq.  (5) 
and  Eq.  (8),  we  see  that  the  interface  shape  can  expressed  in  the  form 


X = 


(9) 


Thus  the  ^-vector  actually  provides  an  explicit  parametric  description  of  the  equilibrium  shape  deter- 
mined by  7(n). 

That  the  argument  of  'y{n)  is  constrained  to  be  a unit  vector  can  be  a complication  when  taking 
variations  of  the  surface  energy;  this  is  one  reason  for  introducing  7(0,  (p).  As  pointed  out  by  Taylor 
et  al.  [9,  10],  a convenient  alternative  is  to  extend  the  domain  of  definition  of  the  function  7(n)  from 
unit  vectors  to  vectors  p of  arbitrary  magnitude  by  setting 


l{p)  = \p\7{p/\p\),  (10) 

which  extends  7 as  a homogeneous  function  of  first  degree.  With  this  extension,  the  function  rj{6,  (p) 
used  in  the  above  definition  (6)  of  the  ^-vector  is  given  by 


r'y{6,(p)  = r'y[n)  = ^{rn), 


(11) 


which  leads  to  the  definition 


(j{p)  = 


d'ljp) 
dpj  ’ 


(12) 


which  is  valid  for  vectors  p of  arbitrary  length.  For  example,  if  the  interface  is  expressed  as  a level  set 
<p{x)  = constant,  then  the  vector  p with  components 


d4> 

dxj 

is  proportional  to  the  normal  vector,  and  the  surface  tension  is  given  by  7 = 7(V^);  Eq.  (12)  then 
defines  |*  = |*(V0). 

Since  j^p)  is  a homogeneous  function  of  first  degree,  it  satisfies  the  identity 

7(?0  = Pj^^  = P-|‘(k).  (14) 

which  results  from  differentiating  the  expression  7(A^  = with  respect  to  A and  setting  A = 1. 

We  note  that  here  and  throughout  this  paper  we  employ  the  Einstein  summation  convention  with 
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the  dummy  indices  ranging  from  one  to  three.  Differentiating  the  expression  7 = p • ^ and  using  the 
definition  Eq.  (12)  to  obtain 

d'r  = C-dp,  (15) 

results  in  the  expression 

p-df=0,  (16) 

Cahn  and  Hoffman  use  the  fundamental  relations  (15)  and  (16)  to  provide  a physical  interpretation  of 
the  {-vector.  They  show  that  — { represents  a force  exerted  by  an  elemental  area  of  the  interface  which 
acts  to  rotate  and  contract  the  element  in  such  a manner  as  to  reduce  its  free-energy.  In  particular, 
the  component  of  { in  the  direction  n acts  to  reduce  the  area  of  the  element  without  rotation  and  the 
component  orthogonal  to  ft  acts  to  rotate  the  element  without  contraction  of  its  area. 

Just  as  the  surface  divergence  of  the  {-vector  replaces  the  term  jK  in  the  Laplace- Young  equation, 
the  isotropic  Gibbs-Thomson  equation 

T=Tm-  ^7^  (17) 

assumes  the  form 

T = Tm-^^s-(  (18) 

in  the  anisotropic  case;  the  main  goal  of  the  present  paper  is  to  estabhsh  that  this  equation  is  recovered 
in  an  appropriate  sharp-interface  limit  of  a wide  class  of  anisotropic  phase-field  models. 

We  note  that  Hoffman  and  Cahn  were  able  to  relate  the  equilibrium  shape  given  by  Eq.  (9)  to  the 
Gibbs- Wulff  construction  that  also  gives  the  shape  of  equilibrium  interface.  In  addition,  Cahn  and 
Hoffman  extended  this  theory,  in  a very  natural  manner,  to  the  situation  where  the  anisotropy  of  the 
surface  energy  is  so  severe  that  the  equilibrium  interface  suffers  missing  orientations  or  facets. 

3.  (^-vector  Formulation  of  Phase-Field  Models 

In  this  section  we  show  how  the  {-vector  arises  in  a natural  way  in  a broad  class  of  anisotropic  phase- 
field  models.  Many  phase-field  models  for  the  solidification  of  a pure  material  consist  of  a thermal 
diffusion  equation  coupled  to  a form  of  modified  Cahn- Allen  equation  whose  double  well  minima  have 
depths  that  vary  with  temperature.  Specific  models  differ  in  the  form  of  bulk  free  energy  density 
that  is  employed,  and  in  the  particular  convention  used  to  normalise  the  bulk  values  assumed  by  the 
phase-field  function.  In  the  sharp  interface  limit,  the  asymptotic  analysis  of  the  thermal  diffusion 
equation  produces  the  usual  interfacial  jump  conditions  for  the  temperature  field,  consisting  of  the 
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continuity  of  temperature  and  the  balance  of  latent  heat  generation  by  the  jump  in  heat  flux  across 
the  interface.  The  modifled  Cahn-Allen  equation  produces  the  condition  for  thermal  equilibrium, 
including  the  Gibbs-Thomson  effect  and  the  effect  of  interface  kinetics.  For  our  purposes  it  suffices  to 
limit  the  discussion  to  the  treatment  of  the  modified  Cahn-Allen  equation  only,  and  for  simplicity  we 
derive  the  equation  from  a relatively  uncomplicated  form  of  Helmholtz  free  energy  functional.  The 
specific  form  of  free  energy  density  that  is  employed  is  unimportant  to  the  discussion,  because  in  these 
models  the  surface  tension  anisotropy  is  produced  through  the  choice  of  gradient  energy  expression  in 
the  free-energy  functional. 

We  consider  a Helmholtz  free  energy  functional  having  the  form 


T 


(19) 


4a  ' Tm  ^ 

where 

= P{^)  = 4,^3-2<t>),  (20) 

77  is  the  gradient  energy  coefficient  and  1/a  is  proportional  to  the  barrier  height  of  the  double- well 
potential.  Here  the  phase  field  assumes  the  values  ^ = 0 and  (f)  = 1 in  the  solid  and  liquid  phases, 
respectively,  and  the  function  P(^)  varies  monotonically  from  P(0)  = 0 to  P(l)  = 1.  The  gradient 
energy  coefficient  77  is  constant  in  the  isotropic  case.  The  governing  equation  for  (f)  is  then  assumed  to 
have  the  form 

i;GW)+  (21) 

where  tti  is  a mobility  that  is  also  constant  in  the  isotropic  case.  For  T = T m the  equation  admits  a 
steady  one-dimensional  solution 


4>(x)  =21 


( 


\ 


\2>/2[77v^y 


+ 1 


(22) 


that  represents  a stationary  planar  solid/liquid  interface  with  a thickness  that  can  be  characterised 
by  the  width 

6 = 2>/2[77\/^.  (23) 


The  interface  has  a surface  energy  7 given  by 


(24) 


^y/2a  ’ 

which  is  proportional  to  the  product  of  the  interface  thickness  and  the  energy  per  unit  volume  1/a 
associated  with  the  double- well  barrier  height. 


-7- 


PHASE-FIELD  ^VECTOR 


October  17,  199^ 


For  T 7^  Tm  the  equation  admits  a one- dimensional  traveling-wave  solution  [16] 

x-vr 


where  the  velocity  is  proportional  to  {Tm  — T), 

6y/2m[7]^/a\L{TM  - T) 


V = 


Tm 


(25) 


(26) 


The  form  of  this  relation  suggests  that  the  mobility  parameter  m can  be  related  to  the  interface  kinetic 
coefficient  ^ by 


M = 


^\/2m[r]^/a\L 

Tm 


(27) 


a result  which  can  also  be  established  through  a sharp-interface  limit  such  as  considered  below. 
The  free  energy  functional  can  be  rewritten  in  the  form 


f = /v  { 


(28) 


where  u = [T  — Tm)ITm  is  the  temperature  relative  to  the  melting  point  in  units  of  Tm-  The  isotropic 
governing  equation  can  then  be  written  in  the  form 

72[aT]  f 'y\  d<l> 


Tm^ 


(29) 


The  sharp-interface  limit  corresponds  to  the  case  that  the  diffuse  interface  width  is  small  compared 
to  a characteristic  macroscopic  length  scale  R,  which  can  be  taken  to  represent  a typical  radius  of 
interfacial  curvature  in  the  nonplanar  case.  We  consider  the  distinguished  limit  r]y/afR  —*  0 and 
aL  0,  while  maintaining  a Unite  value  for  the  ratio  77/  ■y/a  in  order  that  the  surface  energy  remains 
finite  in  the  limit. 

It  is  convenient  to  measure  length  in  units  of  R and  energy  density  in  units  of  T,  in  which  case  the 
free  energy  functional  has  the  form 


= / |36e72|V(/ip  + ^G{(p)  - P{4>)u\  dV,  (30) 

where  e = aL  plays  the  role  of  a small  parameter  in  the  subsequent  asymptotic  treatment.  We  choose 
a diffusive  time  scale  r consistent  with  the  accompanying  thermal  diffusion  processes,  and  scale  the 
dimensional  kinetic  coefficient  in  units  of  RI{tTm)\  the  specific  choice  of  time  scale  is  unimportant 
to  the  discussion  provided  that  the  resulting  dimensionless  kinetic  coefficient  can  be  assumed  to  be  of 
order  unity.  The  dimensionless  governing  equation  then  has  the  form 

= 72e7^VV  - ^G'(0)  + P'{4>)v^.  (31) 

a ot  4e 
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To  formulate  an  anisotropic  model  in  which  the  surface  energy  depends  on  orientation,  Kobayashi 
[4]  defines  a normal  vector  v to  the  level  curves  of  the  phase  field  at  each  point  in  space  by  setting 

V<p 


1/  = 


|V<A| 


(32) 


and  inserts  the  function  7(V0/|V^|)  in  the  relation  (24)  to  obtain  a gradient  energy  coefficient  rj  that 
depends  on  the  orientation  of  the  level  curves  of  (j).  However,  by  expressing  the  free  energy  functional 
directly  in  terms  of  7 as  in  Eq.  (30),  and  employing  the  homogeneous  extension  of  7,  with 

V(j> 


|V0|7 


m 


) = 7(V0), 


(33) 


the  energy  functional  can  be  expressed  in  the  form 

J^=  |36e[7(V</.)]2  + Lg(^)  - P(4>)u^  dV. 

In  this  form  it  is  straightforward  to  take  the  variation  of  the  functional  to  obtain 


(34) 


(35) 


We  note  the  occurrence  of  the  ^-vector,  with  components  = ^7/^^ifcj  in  the  above  expression.  The 
anisotropic  form  of  the  dimensionless  governing  equation  can  therefore  be  written  in  the  form 


= 72.V  . (71)  - 


(36) 


here  we  have  assumed  that  7(^  and  ^{p)  are  homogeneous  functions  of  first  degree  in  order  to  include 
the  effects  of  both  surface  tension  anisotropy  and  kinetic  anisotropy. 

We  note  that  in  this  notation  the  anisotropy  suggested  by  Caginalp  and  Fife  [13]  can  be  expressed 
in  the  simple  form  7(V<^)  = y/{r}i(f)xiy  + (772<^x2)^j  where  771  and  772  are  constants. 


3.1.  Curvilinear  Coordinates 

In  order  to  compute  the  sharp-interface  limit,  it  is  convenient  to  employ  a curvilinear  coordinate 
system  {u^  that  is  based  on  the  assumed  limiting  shape  of  the  solid/liquid  interface  defined  by 

the  level  curve  (j)  = 1/2,  henceforth  denoted  as  5.  At  any  instant  of  time,  the  solid/liquid  interface  for 
e = 0 is  expressed  parametrically  as  X(tt^,u^),  where  v}  and  v?  provide  a local  coordinate  description 
of  the  surface,  as  discussed  in  more  detail  in  Appendix  A.  The  transformation  from  the  curvilinear 
coordinates  to  physical  space  is  given  by 

i(u^,it^,u^)  = X{v} -f  u^n[u^ jU^),  (37) 
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where  n is  the  normal  to  the  interface  at  the  point  The  interface  therefore  corresponds 

to  the  surface  = 0.  The  transformation  generates  local  contravariant  basis  vectors  ej  = dxjdu^ 
and  an  associated  metric  gjk  — ' ^k-  The  partial  derivatives  of  the  transformation  are  denoted  by 

= dx^ Idu^  and  6*  = du^jdx^]  both  are  expressed  as  functions  of  the  variables 
In  the  curvihnear  coordinates  the  phase  held  is  written  as  = <l>{x^,x^,  x^).  To  simplify 

the  notation,  we  denote  the  respective  gradients  by 

d<f>  ^ . . 

which  transform  as  covariant  vector  components  according  to  the  law 

Pj  = b)pk-  (39) 


In  the  curvilinear  coordinates  the  transformed  surface  energy  7 is  then  obtained  by  inserting  the 
expressions  (39)  into  the  function  7(pi,P2jP3);  since  7 then  has  dependence  not  only  on  pk  but  also  on 
through  the  elements  of  the  transformation  matrix  6*  (u^,  u^,  u^),  the  transformed  surface  energy 
is  written  in  the  form 

l{Pi,P2,P3,'u},u'^,u^)  = 7(P1jP2,P3).  (40) 


Because  the  pj  variables  depend  linearly  on  pk,  J is  seen  to  be  a hrst-order  homogeneous  function  of 
the  variables  pk- 

The  transformed  free-energy  functional  takes  the  form 

^ — J |36e7^  -1-  ^G{$)  — P(^)u|  y/gdu^  du^  du^,  (41) 

where  g denotes  the  determinant  of  the  metric  tensor  and  yfg  dv}  du^  du^  is  the  transformed  volume 
element.  The  variation  of  X is  given  by 


8^ 


We  set 


?*  = 


^7 

dpk 


(42) 


(43) 


each  is  a homogeneous  function  of  the  pk  variables  of  degree  zero.  From  the  definition  (40)  it  is 
seen  that  the  variables  form  the  contravariant  components  of  the  transformed  ^-vector;  i.e.,  they 
satisfy  the  relation 

e = • (44) 
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The  expression 


(45) 


is  the  divergence  in  curvilinear  coordinates  (see  Appendix  A)  of  the  product  of  7 and  the  ^-vector. 


4.  The  Limit  e 0. 

In  the  limit  e — ^ 0,  we  expect  the  formation  of  a thin  interfacial  layer  about  the  surface  z)  = 1/2 

in  which  (p  varies  rapidly  between  zero  and  unity.  We  assume  that  the  surface  (p^x^y^z)  = 1/2  moves 
with  normal  velocity  into  the  liquid  phase,  given  by  ^ = 1. 

In  this  anisotropic  formulation  of  the  phase-held  equation,  the  main  technical  difficulty  in  con- 
ducting an  asymptotic  analysis  lies  in  dealing  with  the  term  V • (7^)  in  the  interfacial  layer.  We  now 
examine  its  expansion  in  the  interfacial  layer. 


4.1.  Expansion  of  V • (7^ 

In  the  limit  e —»■  0 we  anticipate  the  formation  of  an  interfacial  layer  of  thickness  0[e).  We  therefore 
introduce  the  rescaled  coordinate  p,  defined  by  v?  = ep  and  represent  the  phase  field  in  this  region  as 

p = -1-  0{e^).  (46) 


The  inner  expansion  of  the  gradient  of  p is  then  found  to  be 

V?!>=  4^°)S+0(l). 

We  rewrite  the  divergence  as  the  sum  of  three  terms, 

V9  ^ Vb  h ^ ■ 

The  inner  expansion  for  this  expression  can  be  obtained  using  the  intermediate  results 

7 = + Itn)  • + 0(«), 

+ 0{e% 


(47) 

(48) 

(49) 

(50) 

(51) 
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that  are  derived  in  Appendix  B;  here  7(n)  and  |*(n)  correspond  to  the  surface  energy  and  vector 
evaluated  at  the  normal  to  the  leading-order  interface  position,  and  are  independent  of  p.  The  surface 
divergence  operator  also  corresponds  to  the  interface  p = 0.  We  then  have  that 


du^  = + [7(«)]^4y  + 7(n)f(") ' + 0(e), 


and 


{V9i'e)  + • {7(»)f(«)4'”}  + 0(e). 


We  now  go  on  to  develop  the  solution  for  (p  at  each  order  in  the  interfacial  layer. 


4.2.  Leading  Order 

The  leading  order  problem  for  is 


72[7(™)]=^  - iG'(0(°))  = 0, 


where  = 40(0  — 1)(0  — 1/2),  with  boundary  conditions 


0(°)^  1 as  p — »•  -t-cjo,  and  0^°^  — + 0 as  p — >■  — oo. 


Hence  0(°)  is  given  by  the  planar  interface  solution  discussed  above  and  so 


0(°)('1£\u2,p)  = - 


tanh 


(247(n)) 


.247(n) 

The  leading-order  dimensionless  surface  energy  is  7(n),  and 

1 


+ 1 


r 

v/  — oo 


= 1 = 


727  (n)  ’ 


from  which  it  follows  that 


VsT= slin)] 

7(n) 


we  shall  use  this  identity  below,  when  invoking  the  solvability  condition. 


4.3.  First  Order 


The  first  order  problem  is 


72[7(™)1"^^  - iG"(</.'°))0'')  = 


72^t;„0W  - 727(jJ)f(n)  ■ - 72  V5  • {7(»)?(«)0',°>}  - (59) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 
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where  we  have  employed  the  homogeneity  of  7 and  Eq.  (47).  Here  —Vn4>^p^  is  the  leading  order  con- 
tribution obtained  by  transforming  the  time  derivative  to  the  curvilinear  coordinate  system.  The 
homogeneous  adjoint  solution  to  the  linearised  operator  on  the  left-hand  side  is  The  solvability 

condition  for  this  equation  gives 

u = ~72  7(n)f(n)  • J dp  ~ 72  j dp  “ 

We  have 

72  7(«)f(n)  • dp  = 36  T(n)f(")  • Vs  dp  = , (61) 

where  we  have  employed  the  identity  (58)  in  the  last  step.  In  a similar  fashion  we  also  find  that 

72  y_”  Vs  • = Vs  • !*(«)  + . (62) 

and  the  solvability  condition  (60)  therefore  simplifies  to  give 

u ~ - Vs  • f(^)  “ VnlpL{n).  (63) 

This  is  our  principal  result.  It  is  the  so-called  Gibbs-Thomson-Herring  equation  expressed  in  term  of 
the  Cahn-Hoffman  ^‘-vector.  It  relates  the  dimensionless  interface  temperature  u to  the  surface  energy 
and  interface  kinetics. 

5.  Discussion 

The  ^-vector  formulation  originally  developed  by  Hoffman  and  Cahn  was  used  in  the  context  of  a 
sharp  interface  model  where  it  is  a vector  field  with  a domain  restricted  to  the  surface  of  the  interface. 
We  have  shown  above  that  the  same  ^'-vector  arises  naturally  in  phase-field  models  that  describe  a 
diffuse  interface.  In  this  setting  the  |-vectoris  defined  throughout  the  whole  domain  as  the  vector 
field  ^(f^),  i.e.,  the  ^-vector  at  a point  in  space  is  given  by  the  Hoffman-Cahn  ^-vector  evaluated  at  the 
local  normal  i/  to  the  level  curve  of  the  phase  field. 

As  noted  in  section  2 in  the  context  of  a sharp  interface  description, 

1{P)  = P-({P),  (64) 

where  the  normal  vector  to  the  interface  is  given  by  n = p/|p|.  This  result  also  holds  in  the  phase-field 
context  in  which  p various  throughout  the  domain,  and  ^ = p/\p\  is  the  local  normal  to  the  phase-field 
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contours.  Further,  we  note  that  along  a curve  we  have  that 

d'y  _ dj  dpk  ^ dp 
duj  dpk  dw  dw 


(65) 


where  dpidw  is  the  covariant  derivative  of  p along  the  curve.  Differentiating  the  identity  (64)  along 
the  curve  and  using  Eq.  (65)  yields 

di 

p-  - = 0,  (66) 

where  d^jduj  is  the  covariant  derivative  of  ^ along  the  curve.  Eq.  (64)  to  Eq.  (66)  are  the  extension  of 
the  fundamental  relations  of  the  vector  given  by  Eq.  (14)  to  Eq.  (16)  appropriate  to  the  phase-ffeld 
interpretation  of  the  ^-vector  . 

Finally,  we  provide  a brief  heuristic  discussion  of  the  relation  of  the  surface  divergence  of  the  ^ 
vector  to  the  Herring  formula  using  the  notation  developed  here.  To  do  so,  we  consider  a non-umbilic 
elliptic  point  [18]  of  the  surface  where  both  principal  curvatures  are  positive.  The  coordinates  are 
taken  to  lie  along  the  lines  of  principal  curvature,  forming  an  orthogonal  grid  near  the  point.  For 
these  coordinates  we  then  have  ti  - 12  = 0,  and  the  derivatives  of  the  normal  vector  are  given  by  the 
Rodrigues  formula  [18] 


dfi 


dfi 


= JC2t2, 


(67) 


du^  du^ 

where  K-i  and  JC2  s-re  the  principal  curvatures  of  the  surface.  This  implies  that  in  the  surface  the 
Christoffel  symbols  satisfy  F^^^^  = —hnKi,  T^22  — “^22^2?  and  = ^^21  — ^ (see  Eq.  (92)  in 
appendix  A).  We  further  assume  that  the  surface  metric  tensor  reduces  to  the  unit  matrix  at  this 
point,  and  that  its  first  derivatives  vanish  there  (corresponding  to  the  surface  coordinates  satisfying  a 
“local  flatness”  condition  at  this  point,  see  [19]).  Under  these  assumptions,  the  surface  divergence  of 
the  ^-vector  at  this  point  takes  the  form 

di^  , 


For  this  point,  it  also  follows  from  Eq.  (102)  in  Appendix  B that  for  p = n and  = 0 we  have 


^7 


= -F' 


dj 


Im  Q -r 


dpr 


dj 

dv? 


= -F' 


^7 


2Tn  Q ~ 


dpr 


= ^2i\ 


SO  that  the  surface  divergence  can  be  written  in  the  form 


d r 1 ^7  1 


d 


^ du'  1 du^  ] du^{  IC2  du^  } 


Defining  new  coordinates  6i{u^)  and  62{u^)  by 


du^ 


4 = ^1, 


du^ 


= }C, 


(68) 


(69) 


(70) 


(71) 
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gives  that 

V..f  (72, 

This  expression  inserted  into  Eq.  (63)  gives  the  extension  of  the  formula  due  to  Herring  [11]  for  the 
temperature  of  a moving  curved  interface. 

6.  Conclusions 

The  gradient  energy  term  in  an  anisotropic  phase-held  model  of  the  general  type  discussed  by  Kobayashi 
[4]  can  be  formulated  in  terms  of  the  surface  energy  j{n).  The  variation  of  the  resulting  free  energy 
functional  is  easily  computed  if  the  argument  of  the  surface  energy  has  been  extended  to  vectors 
of  arbitrary  magnitude  as  a homogeneous  function  of  hrst  degree  [9,  10].  The  resulting  variational 
equation  for  the  phase  held  contains  the  ^-vector  of  Hoffman  and  Cahn  [7,  8].  A formal  asymptotic 
analysis  of  the  sharp  interface  limit  [5,  17]  recovers  the  Gibbs-Thomson- Herring  equation  [11],  in  this 
case  expressed  in  terms  of  the  surface  divergence  of  the  ^-vector. 

We  have  assumed  that  the  surface  tension  anisotropy  is  mild  enough  that  missing  orientations 
do  not  occur  in  the  Gibbs- Wulff  construction.  The  recognition  that  the  phase-held  equations  can 
be  formulated  in  terms  of  the  f- vector  provides  a natural  setting  in  which  to  investigate  the  way  in 
which  a phase-held  model  can  describe  the  situation  in  which  missing  orientations  and  facets  would 
be  present  in  a sharp-interface  model  as  well  for  the  development  of  weak  solutions  and  associated 
numerical  methods.  These  are  lines  of  research  currently  in  hand. 
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A.  Differential  Geometry  Results 


Here  we  summarise  some  results  from  the  differential  geometry  of  surfaces;  for  more  details,  see,  e.g., 

[20,  21]. 

A.l.  The  Surface  Divergence 

We  consider  a surface  expressed  parametrically  in  the  form 


= x[v}y). 


(73) 


where  and  are  general  curvilinear  coordinates  for  the  surface.  Contravariant  basis  vectors  ti 
and  t2  for  the  tangent  plane  to  the  surface  at  any  point  are  given  by  the  partial  derivatives 


dX  , 0,  dX 


The  element  of  arclength  ds  along  the  surface  is  defined  by 


(74) 


ds^  = hiidu^  du^  + 2hi2du^  du^  h22du^  du^, 


1 J..2 


.2  J„.2 


where 


hii  — , 


hl2  = fl21  = tl  ' ^2j 


^22  — ^2  ‘ ^2 


(75) 


(76) 


form  the  covariant  components  of  the  surface  metric  tensor  whose  determinant  is  h = hii/122  — (^12)^  = 
\ti  X t2|^-  The  normal  vector  to  the  surface  can  be  written  in  the  form 


n{u^,u^)  = —r=ti  X ^2- 
y/h 

The  covariant  basis  vectors  form  a reciprocal  basis  set  and  are  given  by 

X 11 

X t2  = - (^22^1  - hi2r2)  , = -^n  X ti  = - [-hi2ti  + hufa)  • 

The  operator  V5  can  then  be  formally  defined  by 


(77) 


(78) 


V5  = i 


•1  9 


+ « 


dv}  du^  ’ 


(79) 


so  a scalar  function  f{y}^vr)  defined  on  the  surface  has  a surface  gradient  given  by 

T-7  X ^1 


(80) 
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which  is  seen  to  be  a covariant  representation  of  the  vector  V 5/.  The  surface  divergence  of  the  position 
vector  is  seen  to  be 


Vs 'X  = t 


_ ri 


dx 


dX 


= P - fi  + f ^ •t2  = 2. 


(81) 


du^  du^ 

More  generally,  it  can  be  shown  [20,  21]  that  the  surface  divergence  of  a vector  A{u^,u^)  written  in 


contravariant  form 


can  be  expressed  as 


A{y}ju'^)  = A^{u^  ,u^)t2  + A'^[u^  ,u^)n, 


2/,.l  „.2' 


„,2' 


(82) 


(83) 


where  JC  is  the  mean  curvature  of  the  interface.  In  particular,  taking  A = n{u^,u^),  this  expression 
shows  that  K is  given  by  the  surface  divergence  of  the  normal  vector. 

A. 2.  Surface-Fitted  Coordinates 

We  next  extend  the  coordinate  system  given  above  to  provide  coordinates  in  a neighbourhood  of  the 
surface.  To  do  this,  we  introduce  a third  coordinate  in  the  direction  of  the  normal  vector,  and  express 
a general  position  vector  x as 


= X{u^ -f-  u^n{y} 

This  generates  the  contravariant  basis  vectors  (see,  e.g.,  [21])  given  by 


(84) 


^ -.2  ..3> 


^ „.2 


dx 

dX 

3 

0 dn 

(85) 

du^ 

du^ 

Avr  = 

dv> 

+ 

dv?-' 

dx 

dX 

3 dn 

0\ 

3 

” du^ 

du^ 

4-  vr—  = 
dv? 

+ 

vr , 

dv?' 

ez[v}, 

y)  = 

-*(  1 

y). 

(86) 

The  y-th  element  of  the  vector  is  denoted  by  u^),  where  = dx^ fdu^.  The  covariant 

components  of  the  associated  metric  tensor  are  given  by 


9jk  — 

and  we  note  that  for  this  particular  coordinate  system  that  • €3  = §2  • ^3  = 0 and  ^3  • ^3  = 1. 
The  differential  operator  V can  be  written  in  the  form 


(87) 


V = e 


* 9 


du^  ’ 


(88) 
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where  here  are  the  covariant  basis  vectors  related  to  the  covariant  basis  vectors  by 

02  X 03  2 X ^3  3 X ^2 


e2  = 


0^  = 


(01-02X03)’  ''  (€2-01X03)’  ''  (€3-01X02)' 

We  denote  the  determinant  of  the  metric  tensor  by  g,  which  also  satisfies  (oi  • 02  X 03)  = ^/g.  For  the 
coordinate  system  considered  here,  we  have  0^  = 03  = n[u^,u^).  The  gradient  of  a scalar  function 

then  has  the  form 


Vcj)  = e 


k d(t> 
du^ 


(90) 


In  applying  V to  a vector  field,  one  must  differentiate  the  basis  vectors  as  well  as  the  coefficients 
The  covariant  derivatives  of  the  covariant  basis  vectors  satisfy 

de^ 


du^ 


+ rUe = 0, 


where  F-^  , are  Christoffel  symbols  defined  by 


= -e-  ■ ^ = ^ 


For  the  coordinate  system  considered  here,  we  have 


^ dfi 

r 3/c  = -ra  • — = 0. 


du^ 


For  a vector  field  A expressed  in  contravariant  form, 

A{y}  = A\v}  ^u^)ejj 

it  can  be  shown  [21]  that  its  divergence  is 

y/g 


(91) 


(92) 


(93) 


(94) 


(95) 


Note  that  for  = 0 we  have  ^1  = ti,  ^2  = ^2?  g = h. 

Comparing  the  expression  (95)  with  the  result  of  applying  (88)  to  the  vector  A = n{u^,u^),  for 
which  A^  = A^  = 0 and  A^  = 1,  gives 


1 dy/g 

y/g  du^ 


= e 


X3=0 


dn 


+ e 


dn 

du^ 


= t 


x3=0 


dn 


At 


dn 

du^ 


= Vs-ft=K:, 


(96) 


where  K is  the  mean  curvature  of  the  surface  = 0 as  discussed  above.  We  therefore  have  that 

_ 1 ^ d[y/hA^{u^ .,,1  o _ dA^ 


\/-A 


du^ 


+ A^(u^,u^,0)  K -}- 


du^ 


x3=0 


= V5  • A(u^,'u^,0)  + 


dA^ 


du^ 


(97) 


t3=0 
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B.  Some  Asymptotic  Results 


In  this  appendix  we  give  derive  asymptotic  representations  of  7 and  f in  the  interfacial  layer. 
The  inner  expansion  of  7 = is  obtained  by  writing 


7 = 1 


Id^  I 2 


du^  ’ du'^  ’ e dp 


= -7  e 


1 ^ d^  d^  I 2 


e ^ \ du^  ’ du^  ’ dp 


,u  ,u  ,ep  , 


(98) 


where  we  have  used  the  fact  that  7 is  first-degree  homogeneous  in  the  partial  derivatives  of  <f).  Inserting 
the  expansion  ^ -f  O(e^)  gives 


7 = - I7  (0, 0,  u^,u^,0)  + e [<^^?7pi  + ^i°2^7p2  + + ^^3]  + ^(e^)}  , (99) 

where  the  partial  derivatives  of  7 are  all  evaluated  with  the  argument  list  (0,0,  From 

homogeneity  we  have 


7 (0,0, 


(100) 


where  the  latter  equality  follows  from  the  observation  that  pi  = p2  = 0 and  ps  = 1 are  the  components 
of  the  normal  vector  n{u^,u^).  It  also  follows  that 


Tpj  (0,0,^(°),u\u^,0)  = 7(n). 


(101) 


A calculation  shows  that 


^ ^ ^ ..1 

duL 


dj  db^j  ^ 
dpj  dvT 


= ^'rnr-^Pk=-T^mrrPk^ 


Therefore 


^7 

du^ 


(0,0,P3,1^^H^0)  = -r^msI^Pa  = 0, 


(102) 


(103) 


since,  as  noted  in  Appendix  A,  T^g  = 0 in  this  coordinate  system.  The  remaining  term  can  be  written 
in  the  form 


^i°^7pi (0> 0)  + 0^?7^(0, 0,  0)  = 1^0)  . 

where  the  leading-order  components  of  the  f- vector  for  j = 1 and  j = 2 are 

= 7pi(0,0,l,‘u^u2,0); 


(104) 


(105) 


here  we  have  used  the  homogeneity  of  7^^. . Together  with  (101),  this  shows  that  are  the  curvilinear 
components  of  the  vector  ^(fi),  so  that  to  leading  order  the  ^-vector  depends  only  on  and  v? . 
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The  expansion  of  7 can  therefore  be  written  in  the  form 

7 = + [f(^)  • + 0{e).  (106) 

The  inner  expansion  of  is  required  through  first  order,  and  is  obtained  by  writing 

= ~ (4„i , , 4>p,  “S  ^p)  ■ (107) 

Expanding  in  e gives 

= 7p,  (0,0,4°^■“^“^0)  + «{0i?7pip,  +^i°'7p2P3  +^1.'’7p3pj  +P7u=p,}  + 0(e^),  (108) 

where  the  second  derivatives  of  7 are  all  evaluated  with  the  argument  list  (0,0,  ,u^ ,0).  As 

shown  above,  the  first  term  is  'y{fi).  Differentiating  the  identity 


% (0, 0,  Ap3, 1t\  14^  0)  =%  (0,0,P3,1^^1i^0) 


(109) 


with  respect  to  A shows  that  7pjp3(0, 0,p3, 0)  = 0.  Finally,  since  7^3(0, 0,p3, tt^, 0)  = 0 as 
shown  above,  we  also  have  that  7u3p3  vanishes  as  well.  Therefore  there  is  no  first-order  correction  for 
and  we  have 

f = 7(^)  + 0(62).  (110) 

We  also  note  that  Eq.  (96)  implies  the  inner  expansion 

= ;c  + 0(e).  (Ill) 

u3=ep 


_LM 

du^ 
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